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ABSTRACT 

We present the results of an investigation of the local escape velocity (Vesc) - line strength 
index relationship for 48 early type galaxies from the SAURON sample, the first such study 
based on a large sample of galaxies with both detailed integral field observations and extensive 
dynamical modelling. Values of Vesc are computed using Multi Gaussian Expansion (MGE) 
photometric fitting and axisymmetric, anisotropic Jeans' dynamical modelling simultaneously 
on HST and ground-based images. We determine line strengths and escape velocities at mul- 
tiple radii within each galaxy, allowing an investigation of the correlation within individual 
galaxies as well as amongst galaxies. We find a tight correlation between Vesc and the line- 
strength indices. For Mgb we find that this correlation exists not only between different galax- 
ies but also inside individual galaxies - it is both a local and global correlation. The Mgb- Vesc 
relation has the form: log(Mgb/4A) = (0.32 zb 0.03) log (Vesc /500km/s) - (0.031 zb 0.007) 
with an rms scatter a = 0.033. The relation within individual galaxies has the same slope and 
offset as the global relation to a good level of agreement, though there is significant intrinsic 
scatter in the local gradients. We transform our line strength index measurements to the single 
stellar population (SSP) equivalent ages (t), metallicity ([Z/H]) and enhancement ([a/Fe]) and 
carry out a principal component analysis of our SSP and Vesc data. We find that in this four- 
dimensional parameter space the galaxies in our sample are to a good approximation confined 
to a plane, given by log (Ksc/SOOkm/s) = 0.85[Z/H] +0.43 log(t/Gyrs) - 0.20. It is surpris- 
ing that a combination of age and metallicity is conserved; this may indicate a 'conspiracy' 
between age and metallicity or a weakness in the SSP models. How the connection between 
stellar populations and the gravitational potential, both locally and globally, is preserved as 
galaxies assemble hierarchically may provide an important constraint on modelling. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: abundances - galaxies: formation 
- galaxies: evolution. 



1 INTRODUCTION 
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In the Hubble classification (Hubble 1936) scheme elliptical and 
lenticular (or SO) galaxies are collectively known as early-type 
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galaxies, and are thought to represent the end-point of many bil- 
Uons of years of evolution. Early-type galaxies exhibit smooth mor- 
phologies, appearing as essentially featureless collections of stars 
on the sky. For many years this simple appearance was thought to 
reflect a straightforward and homogeneous behaviour, both dynam- 
ically and in terms of their stellar populations, across a broad range 
in luminosity and size. More recently observations have shown that 
while in many ways the structure of early-type galaxies is intrinsi- 
cally simple there is a rich diversity in their properties that requires 
a more complex understanding of these objects. Such an under- 
standing will yield important information about the formation and 
evolution of structure in the Universe. 

Many different properties of early-type galaxies are found to 
be well correlated with their luminosities. The earliest correlations 
discovered were those relating global quantities of these galaxies. 
The most luminous galaxies were found to have large half-light 
radii Re (Kormendy 1977), low surface brightnesses within Re, 
(le) and large central velocity dispersions ae (Faber & Jackson 
1976). These correlations can be combined if we plot the mea- 
surements in log cTe, log Re, log(Ie) spacc. In this variable space 
it is found that galaxies are confined to a tight plane, known as 
the Fundamental Plane (Djorgovski & Davis 1987; Dressier et al. 
1987). In this case the value of any one of the variables can be 
calculated once the other two are known - early-type galaxies are a 
two-parameter family. The most luminous galaxies were also found 
to be predominantiy pressure- supported (low V/cr, Bertola & Ca- 
paccioli 1975; Illingworth 1977; Binney 1979), have core surface 
brightness profiles (Kormendy 1987; Lauer et al. 1995; Ferrarese et 
al. 1994, 2006; Faber et al. 1997; Kormendy et al. 2008) and boxy 
isophotes (Bender et al. 1988). In contrast the less luminous galax- 
ies are predominantly rotationally supported (Davies et al. 1983) 
with cuspy surface brightness profiles and discy isophotes. These 
observations hinted at a dichotomy in the early-type population 
(Faber et al. 1997; Kormendy & Bender 1996) but the inclusion of 
two-dimensional kinematics reveals a different, more marked sepa- 
ration into two distinct populations (Emsellem et al. 2007; Cappel- 
lari et al. 2007, hereafter Paper IX and Paper X). Another quantity 
that is found to correlate well with ae is the mass of the galaxy's 
central black hole M. (Ferrarese & Merritt 2000; Gebhardt et al. 
2000), which also correlates with many other galaxy properties in- 
cluding bulge mass Mbuige (Magorrian et al. 1998; Mclure & Dun- 
lop 2002; Marconi & Hunt 2003; Haring & Rix 2004). 

As well as relationships between these global, predominantly 
dynamical quantities there are tight correlations relating stellar 
population parameters. The first known of these was the colour- 
magnitude relation relating the total luminosity of a galaxy to its 
B-V colour (Visvanathan & Sandage 1977). Global colours were 
also found to be well correlated with other galaxy properties, most 
notably central velocity dispersion ae and central absorption line 
strengths for a number of commonly observed absorption indices 
(Bender et al. 1992). That colour and line strength should be tightly 
related is not entirely surprising. The fact that a quantity measur- 
ing a global property of the galaxy (in this case the global colour) 
which is dominated by light from the outer parts of the galaxy 
should be closely related to a quantity measured only in the very 
centre (true for both a and the absorption indices) suggests that the 
behaviour of these properties within a galaxy, as well as between 
different galaxies, must also be confined to a relatively narrow re- 
gion in parameter-space. We shall explore further evidence for this 
idea and it's consequences later in this work. 

There is one further relation most closely related to this work, 
that linking a and the magnesium line strength index (either Mgb 



or Mg2) measured in a central aperture (Burstein et al. 1988). This 
is the tightest and best-studied relation linking a dynamical quan- 
tity a with a quantity depending only on the stellar population, the 
Mg index. Many authors have measured this relation for many hun- 
dreds of early-type galaxies (see e.g. Colless et al. 1999), and while 
the precise gradient and zero-point found for the relation vary from 
author to author the tightness of the relation is common to all stud- 
ies. Any successful model of early- type galaxy formation must ex- 
plain this and the other relations discussed above before it can be 
accepted as accurately describing the formation histories of these 
objects. 

The previously discussed relations are all global ones; some 
local relations have also been studied but only with small samples. 
Franx & Illingworth (1990) found that the local colour in ellipti- 
cal galaxies correlated well with the local escape velocity, Vesc, 
whereas the local colour - local a relation has significantly larger 
scatter. A similar result for the local Mg2-Vesc relation was found 
by Davies, Sadler & Peletier (1993), hereafter DSP93. This depen- 
dence on local parameters, which holds both within a single galaxy 
and between a sample of early-type galaxies begins to suggest that 
a key parameter in the formation and evolution of early -type galax- 
ies is the gravitational potential for which Vesc is a proxy. 

In this work we explore the line strength - Vesc relation for 
the SAURON sample of galaxies (de Zeeuw et al. 2002) for which 
integral-field spectroscopy obtained on the SAURON integral-field 
unit (Bacon et al. 2001) and extensive photometry are available. In 
Section 2 we give details of the SAURON sample, the observations 
and the data reduction process for the photometry and the spec- 
troscopy. In Section 3 we discuss the dynamical modelling of the 
sample from which we derive the potential $ and the escape veloc- 
ity Vesc- In Section 4 we present the resulting line strength- Vesc re- 
lations and translate these to the physical properties age, metallicity 
and alpha enhancement. In Section 5 we consider the implications 
of our results in the context of galaxy formation scenarios. Finally, 
our conclusions are summarised in Section 6. 



2 SAMPLE AND DATA 

2.1 Selection 

The sample of galaxies used in this investigation is the SAURON 
sample of 48 early type galaxies (de Zeeuw et al. 2002), divided 
equally between E and SO morphologies (where the classification 
is taken from de Vaucouleurs et al. 1991). This sample is repre- 
sentative of nearby bright early-type galaxies (cz ^ 3000 kms~^; 
Mb ^ — 18 mag) and is fully described in de Zeeuw et al. (2002). 
The sample consists of an equal number of 'cluster' and 'field' ob- 
jects (where cluster objects are defined as those belonging to the 
Virgo cluster, the Coma I cloud and the Leo group) uniformly cov- 
ering the ellipticity-Ms plane. 

2.2 Data 

The photometric data consists of ground-based photometry ob- 
tained in the F555W filter on the 1.3-m McGraw-Hill Telescope 
at the MDM observatory on Kitt Peak (Falcon-Barroso et al. in 
preparation), supplemented by HST observations where available 
(see Table 1 for the complete list.) A relatively large field of view of 
17.1 X 17.1 arcmin^ was used for the MDM observations in order 
to provide accurate sky-subtraction from the images. The space- 
based observations consist primarily of HSTAVFPC2 imaging or 
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Table 1. Properties of the 48 E and SO galaxies from the SAURON sample used in this paper 



Galaxy 


Type 


Re 


Dist 


Rotator 


i 


(M/L)x 


Band 


Ml 


X-I 


(M/L)/ 


HST 


Quality 


d(log Mgb) 


Name 




(arcsec) 


(Mpc) 




n 




X 


(mag) 


(mag) 




Imaging 


of fit 


d(log Vesc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


NGC 474 


S0"(s) 


29 


32.0 


F 


37 


2.86 




-21.94 




2.86 


814W 


3 


0.66 


NGC 524 


S0+(rs) 


51 


23.3 


F 


20 


5.36 


I 


-22.99 


_ 


5.36 


814W 


1 


0.19 


NGC 821 


E6? 


39 


23.4 


F 


79 


3.58 


I 


-22.32 


_ 


3.58 


814W 


1 


0.43 


NGC 1023 


SBO-(rs) 


48 


11.1 


F 


73 


2.90 


I 


-21.97 


_ 


2.90 


814W 


2 


0.09 


NGC 2549 


S00(r)sp 


20 


12.3 


F 


90 


4.84 




-20.49 


0.65* 


3.64 


702W 


1 


0.41 


NGC 2685 


(R)SB0+pec 


20 


15.0 


F 


76 


1.74 


I 


-20.65 


_ 


1.74 


814W 


2 


0.47 


NGC 2695 


SABOO(s) 


21 


31.5 


F 


48 


5.64 




-21.63 


1.18 


3.80 


_ 


1 


0.30 


NGC 2699 


E : 


14 


26.2 


F 


46 


3.21 




-20.81 


0.61 


2.48 


702W 


2 


0.51 


NGC 2768 


E6 : 


71 


21.8 


F 


90 


5.32 


I 


-22.73 


_ 


5.32 


814W 


1 


0.39 


NGC 2974 


E4 


24 


20.9 


F 


56 


4.79 


I 


-21.94 


_ 


4.79 


814W 


1 


0.39 


NGC 3032 


SABOO(r) 


17 


21.4 


F 


38 


1.99 


I 


-20.41 


_ 


1.99 


814W 


1 


-1.46 


NGC 3156 


SO : 


25 


21.8 


F 


67 


1.46 


I 


-20.43 


_ 


1.46 


814W 


1 


-0.07 


NGC 3377 


E5-6 


38 


10.9 


F 


90 


2.31 


I 


-21.19 


_ 


2.31 


814W 


1 


0.74 


NGC 3379 


El 


42 


10.3 


F 


68 


3.43 


I 


-22.11 


_ 


3.43 


814W 


1 


0.32 


NGC 3384 


SBO-(s) : 


27 


11.3 


F 


66 


1.89 


I 


-21.45 


_ 


1.89 


814W 


3 


0.31 


NGC 3414 


SOpec 


33 


24.5 


S 


60 


4.23 


I 


-22.23 


_ 


4.23 


814W 


2 


0.75 


NGC 3489 


SABO + (rs) 


19 


11.8 


F 


60 


0.99 


I 


-20.99 


_ 


0.99 


814W 


2 


0.25 


NGC 3608 


E2 


41 


22.3 


S 


60 


3.73 


I 


-22.22 


_ 


3.73 


814W 


1 


0.47 


NGC 4150 


SOO(r)? 


15 


13.4 


F 


51 


1.43 


I 


-19.91 


_ 


1.43 


814W 


1 


-0.37 


NGC 4262 


SBO-(s) 


10 


15.4 


F 


26 


8.84 




-20.51 


2.24 


4.08 


ACS/475W 




0.52 


NGC 4270 


SO 


18 


33.1 


F 


70 


4.01 


V 


-21.27 


1.06 


3.01 


606W 


1 


0.02 


NGC 4278 


El -2 


32 


15.6 


F 


40 


4.86 


I 


-22.08 


_ 


4.86 


814W 


1 


0.47 


NGC 4374 


El 


71 


18.5 


S 


60 


4.08 


I 


-23.46 


_ 


4.08 


814W 


1 


0.34 


NGC 4382 


SO+(s)pec 


67 


17.9 


F 


78 


2.58 


I 


-23.39 


_ 


2.58 


814W 




-0.29 


NGC 4387 


E 


17 


17.9 


F 


65 


5.23 




-20.19 


2.22 


2.46 


ACS/475W 


1 


0.21 


NGC 4458 


EO - 1 


27 


16.4 


S 


76 


2.43 


I 


-20.42 


_ 


2.43 


814W 


1 


0.44 


NGC 4459 


SO+(r) 


38 


16.1 


F 


46 


2.86 


I 


-21.86 


_ 


2.86 


814W 


1 


0.31 


NGC 4473 


E5 


27 


15.3 


F 


73 


3.38 


I 


-21.87 


_ 


3.38 


814W 


1 


0.42 


NGC 4477 


SBO(s) :? 


47 


16.5 


F 


26 


6.66 


V 


-21.66 


1.28 


4.09 


606W 




0.08 


NGC 4486 


EO - 1+pec 


105 


17.2 


S 


60 


4.90 


I 


-23.43 


_ 


4.90 


814W 


1 


0.50 


NGC 4526 


SABOO(s) 


40 


16.4 


F 


78 


3.54 


I 


-22.54 


_ 


3.54 


814W 


1 


0.09 


NGC 4546 


SBO-(s) : 


22 


13.7 


F 


70 


6.12 


V 


-21.17 


1.15 


4.23 


606W 


1 


0.46 


NGC 4550 


SBQO : sp 


14 


15.5 


S 


82 


3.38 


I 


-20.50 


_ 


3.38 


814W 


1 


0.41 


NGC 4552 


EO - 1 


32 


15.8 


S 


60 


4.33 


I 


-22.45 


_ 


4.33 


814W 


1 


0.28 


NGC 4564 


E 


21 


15.8 


F 


74 


4.29 




-21.03 


0.65* 


3.22 


702W 


1 


0.48 


NGC 4570 


SOsp 


14 


17.1 


F 


90 


3.39 


I 


-21.42 


- 


3.39 


814W 


1 


0.40 


NGC 4621 


E5 


46 


14.9 


F 


90 


3.95 




-22.80 




3.95 


814W 




0.42 


NGC 4660 


E 


11 


15.0 


F 


68 


3.16 




-20.20 




3.16 


814W 




0.56 


NGC 5198 


El - 2 : 


25 


37.4 


S 


60 


6.48 




-22.04 


0.65* 


4.87 


702W 




0.41 


NGC 5308 


SO-sp 


10 


34.1 


F 


90 


3.56 




-22.13 




3.56 


814W 




0.28 


NGC 5813 


El -2 


52 


31.3 


S 


60 


4.23 




-23.28 




4.23 


814W 




0.19 


NGC 5831 


E3 


35 


26.4 


S 


60 


3.94 




-21.85 


0.58 


3.16 


702W 




0.64 


NGC 5838 


so- 


23 


19.8 


F 


70 


5.22 




-21.73 




5.22 


814W 




0.26 


NGC 5845 


E : 


4.6 


25.2 


F 


75 


3.04 




-20.72 




3.04 


814W 




0.39 


NGC 5846 


EO - 1 


81 


24.2 


S 


60 


5.09 




-23.34 




5.09 


814W 




0.30 


NGC 5982 


E3 


27 


46.4 


S 


60 


3.51 




-23.18 




3.51 


814W 




0.30 


NGC 7332 


SO pec sp 


11 


22.4 


F 


83 


2.78 


V 


-21.44 


1.11 


2.00 


WF1/555W 




0.16 


NGC 7457 


SO-(rs)? 


65 


12.9 


F 


64 


1.89 


I 


-20.73 




1.89 


814W 




0.07 



Notes: Column (1): NGC number. Column (2): Morphological type from RC3. Column (3): Effective (half-hght) radius Re measured in the I band (see Paper 
IV). Column (4): Distances were taken from preferentially from Mei et al. (2007) or Tonry et al. (2001). Virgo galaxies without other distance determinations 
were assigned the mean Virgo distance of 16.5 Mpc from Mei et al. (2007). Distances for other galaxies were taken from Paturel et al. (2003). Column 
(5): Galaxy classification from Paper IX: F = fast rotator (Xji > 0.1), S = slow rotator (Xji ^ 0.1). Column (6): The best-fitting inclination determined 
from axisymmetric Jeans Anisotropic MGE (JAM) modelling. For slow rotators an inclination of 60° is assumed (except NGCs 4458 and 4550 which have 
independent determinations of their inclination). Column (7): The M/L of the best-fitting JAM model, in the given band. Column (8): Photometric band the M/L 
was determined in. Column (9): Total magnitude determined from the MGE models and converted to /-band using colours from the literature. Column (10): 
Galaxy colour, where X is the band of the HST imaging. Colours were taken preferentially from Tonry et al. (2001) or from Prugniel & Heraudeau (1998). For 
those galaxies marked with a * no colour was available and an average early-type colour from Prugniel & Heraudeau (1998) was used. Column (11): The M/L 
of the best-fitting JAM model, converted to /-band. Column (12): The instrument and filter on the HST from which the photometry is taken. Unless otherwise 
stated data was taken using HST/WFPC2. Column (13): The quality of fit of the MGE model. For galaxies ranked 1 a good fit was achieved. For galaxies 
ranked 2 the fit achieved was still good apart from minor discrepancies or dust absorption. For those ranked 3 significant discrepancies had to be taken into 
account in the fitting process and the MGE model does not closely follow the isophotes (mostly due to bars) but still reproduces reasonable second-moment 
velocity maps. Column (14): Mgb - Vesc gradient determined as described in Section 4.2. 
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Figure 1. Photometry of NGC4570 with the MGE model contours over- 
plotted. From top to bottom the figures are: HST/PCl field, HSTAVFPC2 
mosaic image, wide-field MDM image. The MDM image is oriented with 
north to the top and east to the left. The WFPC2 mosaic field is shown 
on the MDM image. As can be seen the MGE contours closely follow the 
isophotes at all radii. 



imaging from ACS or WPFC where WFPC2 data was not avail- 
able. 

The HST data were used as reference for the photometric cal- 
ibration and the MDM images were rescaled to the same level. The 
method of photometric calibration is described in detail in Cappel- 
lari et al. (2006, hereafter Paper IV), but in summary we measured 
logarithmically sampled photometric profiles using circular aper- 
tures for each image after masking bright stars or galaxies. We do 
not expect or observe strong colour gradients between the F555W, 
F814W and intermediate filters, allowing us to match the MDM 
and HST photometry. The photometric profiles were then fitted by 
minimising the relative error between the two profiles in the region 
of overlap. The HST and MDM images were then merged to form 
a single photometric profile for each object. 

The spectroscopic information was obtained using the 
SAURON integral-field unit on the 4.2-m William Herschel Tele- 
scope at the Roque de los Muchachos observatory on La Palma. For 
details of the instrument and the data reduction pipeline see Bacon 
et al. (2001). The SAURON field of view covers objects out to typ- 
ically 1 Re and at least 0.5 Re. The data reduction steps include 
bias and dark subtraction, extraction of the spectra using a fitted 
mask model, wavelength calibration, flat fielding, cosmic-ray re- 
moval, sky subtraction and flux calibration. The flux calibration is 
described in detail in Kuntschner et al. (2006, hereafter Paper VI). 
The stellar absorption lines are also properly corrected for nebular 
emission. The SAURON wavelength range allows us to measure 
four Lick indices (Trager et al. 1998): H/3, Fe5015, Fe5270 and 
Mgb (see Worthey et al. 1994, for a full definition of these indices). 
In this work we consider only H/3, Fe5015 and Mgb as the Fe5270 
line lies at the edge of SAURON 's spectral range and so has incom- 
plete spatial coverage in some objects. The measurement of the line 
strength indices from the final data cubes is described in Paper VI 
where line strength maps for the whole sample are presented. The 
stellar kinematics^ we use in this paper is the same that was used 
in Paper IV which was presented in Emsellem et at. (2004). This 
makes our M/L values directly comparable with those of paper 
IV, when scaled to the same distances. 



3 DYNAMICAL MODELLING 

3.1 Multi- Gaussian Expansion mass models 

Photometric models for all 48 galaxies in the sample were con- 
structed using the Multi-Gaussian Expansion (MGE) parametriza- 
tion of Emsellem, Monnet & Bacon (1994). The observed sur- 
face brightness profile is described in terms of the sum of two- 
dimensional Gaussians, which allows the reproduction of ellipticity 
variations and strongly non-elliptical isophotes. The MGE fitting 
method of Cappellari (2002) was used to facilitate fitting of a large 
sample of galaxies. The MGE models were constrained to have 
constant position angle (PA) to enable axisymmetric Jeans mod- 
elling to be used in determining the underlying potential. 

Twenty four of the MGE models used in this investigation 
were discussed in Paper IV and will not be discussed further here. 
The remaining 24 early-type galaxies of the SAURON sample 
are those for which either accurate Surface Brightness Fluctua- 
tion (SBF) distances were unavailable, WFPC2/F814W data was 
unavailable or the objects show strong non-axisymmetric features. 



Available from http://www.strw.leidenuniv.nl/sauron/ 



While accurate distances are required to estimate mass to light ra- 
tios they are not required for this investigation and by relaxing the 
requirement for WFPC2/F814W data we can now model the entire 
SAURON early-type sample. MGEs for those galaxies not already 
presented in Paper IV are listed in the appendix. Note that several 
of the galaxies already modeled in Paper IV are triaxial objects. 

The MGE models were fitted simultaneously to the wide- 
field MDM images and the higher resolution HST images. Where 
WFPC2 imaging was available the models were fitted simultane- 
ously to the ground based, lower resolution mosaic and higher res- 
olution WFPC2/PC1 images. The MGE fits were performed by 
keeping the position angle (PA) of the Gaussians constant and also 
taking the point spread function (PSF) into account. PSFs were cal- 
culated using Tiny Tim (Krist 1993) and modeled using the above 
MGE fitting method. The PSFs used are presented in Table Bl. 
The resulting analytically deconvolved MGE models are all cor- 
rected for galactic extinction following Schlegel, Finkbeiner & 
Davis (1998), as given by the NASA/IPAC Extragalactic Database 
(NED). They are then converted to a surface density in solar units 
in the Johnson-Cousins magnitude system using the calibration rel- 
evant to each instrument (WFPCl - Harris et al. (1991); WFPC2 - 
Dolphin (2000); ACS - Sirriani et al. (2005).) Absolute magnitudes 
for the Sun (Mb = 5.48,Mv = 4.83, Mr = 4.42, Mi = 4.08,) 
are taken from Table 2.1 of Binney & Merrifield (1998). The val- 
ues of the MGE parameterizations are presented in the Appendix 
in Table Dl and D2. 

Because the SAURON field-of-view is relatively small when 
compared to our imaging we are principally interested in fitting the 
central regions of each galaxy, while the MDM imaging is used 
to provide additional constraints on the MGE model. We do not 
attempt to accurately model structure such as shells or isophotal 
twists in the outer parts of these galaxies but we do seek to repro- 
duce the overall shape of the object. As discussed in Cappellari 
(2002) the models were regularized by requiring the axial ratio of 
the flattest Gaussian to be as round as possible while still repro- 
ducing the observations. This is important so as not to artificially 
constrain the possible inclinations of the models and to reproduce 
realistic densities. We also masked dusty regions in the images; in 
the small number of galaxies in our sample that exhibit dust, the 
dust is only visible in one half of the galaxy image and so does not 
reduce the quality of our MGE fits. 

The quality of the resulting models with respect to the pho- 
tometry was visually inspected for all galaxies to ensure a reason- 
able fit had been achieved. We also compared the resulting kine- 
matics (see Section 3.2) to the SAURON kinematics presented in 
Emsellem et al. (2004) and adapted the MGEs (within the rules out- 
lined below) where necessary to obtain a match to the SAURON 
observations. The models were refined until a satisfactory qualita- 
tive fit was achieved for all galaxies. For those galaxies with regular 
photometry this was easily achieved. An example of the model and 
data photometry for such a galaxy, NGC4570 is given in Fig. 1 . 

3.1.1 Bars, twists and non-axisymmetry 

For those galaxies with non-axisymmetric features such as bars or 
isophotal twists achieving a good match proved more challenging. 
This is because by simply fitting the two-dimensional isophotes we 
infer a three-dimensional distribution of matter than does not re- 
flect that in the real object. In these cases a simple prescription was 
followed to produce our MGE models. 

We note that bars are always associated with a discy structure 
which is to first order axisymmetric and of which the bar repre- 
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Figure 2. HST/PCl photometry on NGC3384 showing two different MGE 
model contours overplotted. In the upper figure we have allowed the model 
to follow the isophotes such that the effect of the bar is included in the 
MGE model. In the lower image we have constrained the axial ratios of 
the Gaussians as described in the text in order to reduce the influence of 
the bar on the MGE model. Despite being a superficially poorer fit to the 
photometry the MGE model from the lower figure reproduces a better fit to 
the observed kinematics than that from the upper figure (see Fig. 3). 

sents a perturbation. Because of this we assume that a reasonable 
axisymmetric mass model of a barred galaxy is found when the el- 
lipticity over the barred region is fixed at the global ellipticity. In 
galaxies where there is an obvious bar we constrained the MGE 
model such that the axial ratio of the Gaussians over the barred 
region was consistent with the axial ratios of the inner and outer 
regions where the impact of the bar on the photometry was negligi- 
ble. An example of a model for a barred galaxy is given in Figs. 2 
and 3. With the barred MGE the axisymmetric JAM model fails to 
reproduce the kinematics, we find a significant improvement when 
using our bar-less MGE. Moreover, the recovered inclination of the 
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Figure 3. Bi- symmetrized and linearly interpolated maps of the second velocity 
moment ii'^ = VV^+~^ (left) and velocity field V (right) of NGC3384. The 
top image shows the observed fii^ and V extracted from the SAURON stellar 
kinematics. The centre and bottom images show the predicted and V from 
the two MGE models described in the caption to Fig. 2. The centre image shows 
the predictions for an MGE that follows the isophotes of the bar whereas the 
bottom image shows the predictions for an MGE that attempts to avoid the 
effect of the bar. As can be clearly seen the 'bar-less' MGE is a much better 
fit to the observed and V. The contours shown are the isophotes from the 
reconstructed SAURON image. 

model becomes closer to the one inferred from the disc. This sug- 
gests that our bar- less MGEs is a better approximation to the global 
galaxy structure. 

For galaxies where the photometric and kinematic PA (see ta- 
ble 1, Paper X) differ significantly the kinematic PA was used for 
the MGE models as being more representative of the galaxy over 
the SAURON field of view. Although in this way we do not rep- 
resent the isophotal twists in the photometry the stellar kinematics 
are fitted better leading to a significant improvement in the value 
of the fit. We checked the kinematics produced by the axisymmet- 
ric Jeans modelling and in all cases a reasonable agreement was 
found with the observed kinematics. Further examples of the ob- 
served and modelled first and second velocity moments are shown 
in the Appendix in Fig. CI. 

3.2 Jeans Anisotropic MGE (JAM) axisymmetric modelling 

In order to compute Vesc we constructed Jeans Anisotropic MGE 
(JAM) axisymmetric models^ (Cappellari 2008) of all the galax- 
ies in the sample. For a given inclination i, the MGE surface den- 
sity can be deprojected analytically (Monnet, Bacon & Emsellem 

^ Available from http://www-astro.physics.ox.ac.uk/~mxc/idl/ 



1992) to obtain the intrinsic density z) in the galaxy merid- 
ional plane, still expressed in terms of Gaussians. This deprojection 
is non-unique but represents a reasonable choice as the resulting in- 
trinsic density resembles observed galaxies for all lines of sight. We 
then apply JAM modelling to the resulting deprojected densities to 
determine the underlying potential. 

The method is described fully in Cappellari (2008) but we will 
briefly summarise the key points here. The positions x and veloci- 
ties V of a large system of stars can be described by a distribution 
function /(x, v) which in a steady state must satisfy the collision- 
less Boltzmann equation. In order to make use of this equation fur- 
ther simplifying assumptions must be made. A typical first choice is 
to assume axial symmetry, which leads to the two Jeans equations 
(Jeans 1922) , but this is not sufficient to specify a unique solution. 
We make the further assumptions that the velocity dispersion ellip- 
soid is aligned with the cylindrical coordinate system (i?, 0) and 
that the anisotropy is constant. We also assume that mass follows 
light, but allow for a constant dark matter fraction. Under these as- 
sumptions the Jeans equations reduce to: 

buv^-iyvl d{hvvi) 

— H ^ — —V (1) 

R OR dR ^ ^ 

dz dz 

where b quantifies the anisotropy, v'j^ — bvl, v is the luminosity 
density and Vi the components of the velocity. T hey p rovide the 
second moments of the line-of-sight velocity v^^, = vf^^, which are 
generally considered to be good approximations to the observed 
quantity Vrms = + cr^. By comparing the observed and Jeans 
modelled second moments we determined the best fitting inclina- 
tion, anisotropy and constant mass-to-light ratio M/L for all 48 
galaxies in our sample. 

Determining the inclination is difficult but we apply several 
independent checks which validate our fitted results. For highly 
flattened objects they must be close to edge on (16 objects). Six 
galaxies have an obvious embedded gas disk or dust lane, inclina- 
tions were estimated assuming these are thin discs. In all cases our 
inclinations were consistent with the independent determination. 
For the remaining 26 galaxies the inclinations are determined based 
purely on the model and may not always be accurate. The accu- 
racy of these inclinations is discussed further in Cappellari (2008). 
Our Vesc is only weakly dependent on the inclination used and so 
this uncertainty does not significantly affect the conclusions of this 
work. As an extreme test of the dependence of our modeled Vesc 
on inclination we artificially set the inclinations for the models of 
all our galaxies to 90° (i.e. edge-on) and re-calculated the Vesc- 
The effect on Vesc was small, only a 5 percent change in the most 
extreme cases.. 

As mentioned above some of the objects in our sample are 
clearly not axisymmetric systems and so the use of axisymmetric 
models requires some justification. The alternative would be to use 
the more general Schwarzschild (1979) models. In Paper IV we 
compared the mass-to-light ratios (M/L) derived from axisymmet- 
ric Jeans and Schwarzschild modelling and find an excellent agree- 
ment between the two (see particularly Fig. 7 from that paper.) Ad- 
ditionally, the slow-rotators, while likely to be triaxial objects are 
also very round (see Paper X, van den Bosch et al. 2008) and so any 
deviations from axisymmetry in their intrinsic shapes are relatively 
small. From this one should not expect major biases in the M/L 
we derive with axisymmetric models. We explicitly tested whether 
this is the case, using the M/L derived via more general triaxial 
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Figure 4. Vesc and Mgb maps for several galaxies with the elliptical annuH used to extract the profiles shown. Here we show only a few of the ellipses to avoid 
cluttering the plot. The numbers to the bottom right of each plot indicate the range of values displayed. As can be seen the modelled Vesc held (or equally the 
potential <l>) closely traces the ellipses but the observed Mgb maps are less regular. 



models by van den Bosch (2008) for eight slow-rotators in com- 
mon with our sample. We found good agreement between our M/L 
determinations with JAM and the triaxial models. A more detailed 
comparison will be presented elsewhere. As for the M/L, we expect 
Vesc to be only weakly affected by the assumption of axisymmetry 
and the use of Jeans models. This can be understood by noting that 
while there are many more possible orbits in a triaxial system than 
in an axisymmetric one this does not affect the potential and hence 
Vesc- It is the distribution of the mass, not the structure of the or- 
bits, that affects and this is largely unchanged between triaxial 
and axisymmetric systems, apart from a small geometric factor. 

3.3 Extraction of line strength indices and Vesc 

In order to study the index-Vesc relations we must extract the in- 
trinsic line strength indices from the SAURON maps and the Vesc 
from our JAM models in a consistent fashion. The potential $ is 
calculated as in Emsellem, Monnet & Bacon (1994) and the Vesc is 
simply related to this by the expression: 

Vesc = ^/2\^{R,z)\ (3) 

The observed indices on the sky plane are the luminosity- 
weighted average of the local values in the galaxy along the line- 
of-sight. We make the quite general assumption that the indices are 
related to Vesc by a simple power-law relationship of the form: 

Index oc Vic (4) 

This leads to: 




p{R, z)Vle(i?, z)dz' 




p{R,z)\2^{R,z)P^^dz' (5) 



With this assumption it is possible to extract the luminosity- 
weighted average, Vesc,p of the local Vesc along the line-of- sight. 
In practice the Vesc values depend only weakly on the parameter 7 
and our conclusions hold for any reasonable choice of the parame- 
ter. 

To form profiles we sum the local line-of- sight integrated val- 
ues over elliptical annuli aligned with the kinematic major axis of 
the galaxy and evenly spaced in radii over the entire SAURON field, 
where the ellipticity used is the global ellipticity as given in table 
1 of Paper X. The noise in each profile is minimised by choos- 
ing the photometric ellipticity e for the profile extraction ellipses. 
Several examples of the Mgb and Vesc maps with the elliptical an- 
nuli used to extract the profiles overplotted are shown in Fig. 4. 
Before this was done the individual SAURON line strength maps 
were inspected for irregular bins. These occasionally occur in the 
outer parts of the SAURON field due to the continuum effects de- 
scribed in Paper VI. Masks were constructed for several of the most 
affected maps. Only the outer few elliptical annuli are affected by 
this issue and the use of un-masked maps does not significantly ef- 
fect the Index-Vesc profiles. Bright stars and obvious dust features 
were also masked on the line strength maps. The error in the line 
strength for each data point is the rms sum of the measurement error 
from Paper VI and the rms scatter within an annulus. We adopted 
an error of 5 per cent in Vesc (see Paper IV). 



4 RESULTS 

In this section we present the Index-Vesc relations determined us- 
ing the method described above. The profiles are presented in Fig. 
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Figure 5. The line strength index versus Vesc relations. All indices are measured in A of equivalent width. In each figure all the points for an individual galaxy 
are represented by a single combination of colour and symbol. In all three cases the correlation is remarkably tight. The black line in each panel is a fit to 
the central Re/8 aperture values, excluding outlying galaxies. The profiles with their NGC numbers shown are those galaxies showing evidence of recent star 
formation. Typical error bars are shown in each figure. 
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Table 2. Linear fit parameters for the Index- Vesc relations from Re/8 circular aperture values. 



Sample 


Index 


a 


b 


Standard 










Deviation 




Mgb 


-0.30±0.09 


0.32±0.03 


0.033 


All (44) 


Fe5015 


0.51 ±0.09 


0.07±0.03 


0.033 




H/3 


1.18ib0.15 


-0.33±0.05 


0.049 




Mgb 


-0.27±0.10 


0.31±0.04 


0.030 


'Clean' (34) 


Fe5015 


0.45±0.10 


0.09±0.04 


0.031 




H/3 


1.04=b0.16 


-0.28±0.06 


0.046 




Mgb 


-0.31±0.13 


0.33±0.04 


0.034 


Fast Rotator (32) 


Fe5015 


0.60±0.11 


0.04±0.04 


0.039 




H/3 


1.27ib0.18 


-0.36±0.06 


0.047 




Mgb 


-0.28±0.16 


0.31±0.05 


0.030 


Slow Rotators (12) 


Fe5015 


0.24±0.13 


0.15±0.05 


0.035 




H/3 


0.91 ±0.24 


-0.24±0.08 


0.052 



Notes: We fit a straight line of the form log Index = a-\-b x log Vesc • The linear fit parameters were calculated using a minimisation technique as 

described in the text. 



5. Mgb and Fe5015 show a remarkably tight correlation, with the 
Mgb-Vesc relation having the tightest correlation. The H/3-Vesc 
correlation is less tight. In this and all following linear fits we fitted 
a linear relation to the Re/8 circular aperture values for each of the 
observed correlations by minimizing the parameter using the 
FITEXY routine taken from the IDL^ Astro-Library (Landsman 
1993), which is based on a similar routine by Press et al. (1992) 
and adding quadratically the intrinsic scatter to make x^/^ = 1' 
where ly is the number of degrees of freedom. For a discussion of 
the technique and its merits see Tremaine et al. (2002). Four out- 
lying galaxies were excluded from these fits (see Section 4.2 for a 
description of how outliers were selected). The results of these fits 
are plotted as the straight line on each figure, and the zero-point, 
slope and rms scatter for each relation are shown in Table 2. Mgb 
and Fe5015 rise with increasing Vesc, with Mgb having the steeper 
slope. In contrast H/3 shows the opposite trend, in the sense that the 
areas of deepest potential $ have the weakest H(3 absorption. 

It is known that some early type galaxies may be weakly triax- 
ial (Kormendy & Bender 1996, Paper IX; Paper X) and an axisym- 
metric model may not reliably reproduce the intrinsic kinematics 
and potential. Galaxies where this is the case typically show a large 
kinematic and photometric misalignment; barred galaxies are also 
not fully axisymmetric systems. Some of the scatter observed in 
the Index- Vesc relations may be due to axisymmetric models not 
properly reproducing the intrinsic Vesc- In order to investigate this 
we define a 'clean' sample in which all galaxies that show evidence 
for a non-axisymmetric distribution have been removed (see Table 
1, column (13), only galaxies graded 1 were included in this clean 
sample.) The best-fitting linear fit parameters for the clean sample 
are given in Table 2. There is little change in the gradients between 
the full and axisymmetric samples, though the scatter is slightly re- 
duced by 10 per cent. This suggests that while imperfect fitting 
of $ due to the assumption of axisymmetry accounts for some of 
the scatter observed in the relations it is only a small effect. 

In Emsellem et al. (2007) a classification scheme 
was described for galaxies based on a parameter Xr = 
{R\V\)/ {R^/V^ + cr^) which is related to the angular mo- 
mentum per unit mass of the stars integrated within 1 Re. Within 
this classification galaxies with Xr < 0.1 are described as slow 

^ http://www.ittvis.com/ProductServices/IDL.aspx 



rotators and those with Xr > 0.1 as fast rotators. Emsellem et 
al. (2007) and Cappellari et al. (2007) speculated that slow and 
fast rotators represent two different families of galaxies with 
significantly different formation histories (in terms of interaction 
or merger events, cold gas accretion episodes, secular evolution 
etc.). If this is the case we might expect the stellar populations of 
the two galaxies to have experienced different histories and for 
signatures of this to show up in the line strength- Vesc relations. 
Definite predictions of these differences are beyond the scope of 
this work but it is thought that dry merger processes are more 
important in the formation of slow rotators whereas the role of 
gas is more prominent in fast rotators. Mergers would tend to 
alter Vesc while leaving the stellar population (and therefore 
the line strengths) unchanged whereas gaseous processes will 
alter the line strengths. In order to explore this we separate our 
sample into fast- and slow-rotators and again perform linear fits 
to the two sub-populations; the best fitting parameters are again 
shown in Table 2. For the case of Mgb the fast- and slow-rotators 
follow the same relationship. In the case of H/3 and Fe5015 there 
is some suggestion that fast- and slow-rotators follow different 
relationships, though when only 'clean' galaxies are considered 
this disappears. 

We also explored the dependence on the traditional divi- 
sion into ellipticals/SOs based on their RC3 classifications (though 
Emsellem et al. (2007) and Cappellari et al. (2007) suggest the 
slow/fast rotator classification is more physically relevant) with 24 
of each in our sample and found no significant difference in the 
Index- Vesc relations between the two sub-samples. We also split 
our sample into field galaxies and those belong to a cluster or group 
(again with 24 galaxies lying in each sub-sample) and again found 
no significant differences but it important to remember that the en- 
vironmental classification used here is a simple one. 



4.1 Comparison to previous work 

The number of studies that have looked at the local Mg-cr and Mg- 
Vesc relations is surprisingly small given the well known tightness 
of the global relation. The two main studies in the area are DSP93 
and CaroUo & Danziger (1994). Both have much smaller samples 
than our current work (8 galaxies and 5 galaxies with Vesc respec- 
tively). The DSP93 sample has four galaxies in common with the 
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Figure 6. Comparison of the results of our sample with those galaxies we 
have in common with DSP93. The Mg2 index values of DSP93 and CD94 
were converted to Mgb index values using Equation 6. The open symbols 
are the values from DSP93 and the closed symbols are the results from this 
work. The solid line is the fit to the SAURON sample, the dotted line to the 
DSP93 sample and the dashed line the CD94 sample. 



SAURON sample ( NGC 3379, 4278, 4374 and 4486) for three of 
which DSP93 have Vesc (the exception is NGC 4374) whereas we 
have no galaxies in common with the CD94 sample. 

Both DSP93 and CD94 looked at the Mg2-Vesc relation, 
rather than Mgb-Vesc relation. Mg2 is a broader molecular index 
but is tightly correlated with the Mgb index (J0rgensen 1997). We 
convert their Mg2 index values to Mgb using: 



log Mgb = 1.57 Mg2 + 0.208 



(6) 



DSP93 and CD94 are also based upon long-slit spectroscopy 
rather than integral-field data. In order to fairly compare our data to 
the previous work we re-extract our Vesc, cr and Mgb profiles using 
a rectangular aperture 2.5 arcseconds wide (as used in the DSP93 
observations) aligned with the major axis of the SAURON maps, 
and sampling linearly in distance along the slit from the centre of 
each galaxy. 

For the three galaxies in common with the DSP93 sample 
we find reasonable agreement with our Mgb-Vesc result (see Fig. 
6). While the individual measurements are in broad agreement 
the slope for our sample is significantly different to that found 
by DSP93 and CD94. This is largely because we sample a much 
broader range in Vesc, approximately twice that in DSP93 and 
CD94. 
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Figure 7. Upper panel: Mgb vs. Vesc measured within circular Re/8 aper- 
tures. The solid line is a linear fit to the data. The sohd symbols mark the 
outlying galaxies which have been excluded from the fit. Lower panel: The 
gradients of the Mgb-Vesc relations for each individual galaxy profile, de- 
termined by fitting a straight line to the galaxy profiles in the same way 
as for the global relation. The histogram shows the individual galaxy gra- 
dients. The dotted line is a continuous distribution with the same mean, a 
and total area as the individual gradients. The vertical solid line shows the 
global gradient determined from fitting to Re/8 values only, with the dashed 
lines indicating the 2 cr error. 



4.2 The gradients within galaxies 

The Mgb-Vesc relation is particularly interesting, partly because it 
is the tightest correlation but also because the profiles for individual 
galaxies follow the global relation remarkably closely. This agrees 
with the results found by DSP93 and CD94. To better illustrate this 
important point we performed a linear fit to each of the individ- 
ual galaxy profiles. In Fig. 7 we show the gradients determined by 
a linear fit for each galaxy, along with the global gradient deter- 
mined from a linear fit to the Re/8 circular aperture value for all the 
galaxies. The global gradient is 0.32 zb 0.03. The distribution of the 
individual galaxy gradients has a biweight mean of 0.34 and robust 
cr of 0.2 (see HoagUn, Mosteller & Tukey 1983, for a description of 



robust statistics). The typical error in the individual galaxy gradi- 
ents is 0.04. The mean of the individual gradients is consistent with 
the global gradient within the errors, but the distribution is con- 
siderably broader. The additional observed scatter in the individual 
gradients implies an intrinsic scatter of 0.16. The Fe5015 and H^^ 
relations behaves quite differently. In the case of Fe5015 the local 
gradients are typically steeper than the global gradient. The local 
gradients in Mgb and Fe5015 vs Vesc appear the same, whereas 
the global gradient is significantly flatter in the case of Fe5015. 
Galaxies typically show H/3 to be flat or slightly rising with Vesc, 
whereas (as mentioned above) the global trend is that H/? falls with 
increasing Vesc- This is in line with studies of radial gradients in 
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Figure 8. The residuals of the Mgb-Vesc relation plotted against llf3. There 
is a clear trend with H/3, in the sense that the residuals are larger (in the 
absolute sense) as H/3 increases. This trend is largely driven by those galax- 
ies that stand out in the Index- Vesc relations, shown as red symbols. The 
majority of galaxies cluster around a small region centred on zero residual. 
The straight line is a fit to the red points only. 

early-type galaxies (e.g. Mehlert et al. 2003) which find galaxies 
have typically very uniform ¥1(3 indices and hence characteristic 
ages. 

Taking a 2a cut in Fig. 7 we note that 3 galaxies have signif- 
icantly different gradients from the mean: NGCs 3032, 4150 and 
4382. These three galaxies also stand out in the Mgb-Vesc relation, 
deviating significantly from the mean relation. NGC 3156 also de- 
viates significantly from this relation, and even though its local gra- 
dient lies within our 2 a cut it has the 2nd largest error on it's local 
gradient due to it's U-shaped profile. For this reason we also con- 
sider NGC 3156 as an outlier. These four galaxies are labelled in 
Fig. 5 with their NGC numbers. Three of these galaxies have the 
highest values of H/3 in our sample, indicative of recent star forma- 
tion (see Paper VI). The fourth, NGC4382, is a peculiar galaxy in 
both the kinematic and line strength maps, showing a central dip in 
a and Mgb as well as a disturbed morphology. It also has signif- 
icantly higher H;^ than galaxies of a similar Vesc, associated with 
star formation in a central disc. While these galaxies stand out no- 
ticeably in Mgb and ¥1(3 they lie much closer to the Fe5015 relation, 
with only NGC3156 showing a significant deviation. 

Three of these outliers also have the lowest values of Vesc in 
our sample. It is possible that there is a break in the Mgb-Vesc re- 
lation at these low values of Vesc but we cannot make a definitive 
judgement on this, given the limited number of galaxies in our sam- 
ple that occupy this regime. Either low- Vesc galaxies are simply the 
least massive galaxies, expected to have experienced more recent 
star formation in a downsizing scenario, in which case we might 
expect them to return to the observed relations, or the Mgb-Vesc 
relation breaks down at these low values of Vesc, suggesting dif- 
ferent processes determine the stellar population characteristics of 
galaxies in this regime. A sample with more galaxies in this regime 
would be required to decide between these two hypotheses. 

4.3 Accounting for H/3-strong galaxies 

We mentioned above that those galaxies that deviate significantly 
from the Mgb-Vesc relation also have high values of H/3. We 



Figure 9. The Index- Vesc relation corrected using the relationship between 
the Mgb residuals and ¥1(3 shown in Fig. 8. The remaining scatter in this 
relation is consistent with the measurement errors. 



quantify this in Fig. 8 by plotting the Mgb residuals, AMgb = 
Mgb(observed) - Mgb(fitted), against H/3. While most galaxies 
cluster in a large clump centred on AMgb = there is a signif- 
icant tail of points with large residuals which appear to correlate 
with H/3. A linear fit to only those galaxies with high values of H/3 
gives the relationship: 

A log Mgb = (-0.36 it 0.05) log R(3 + (0.07 ± 0.03) (7) 

We can use this relationship to modify our Mgb-Vesc relation for 
H;^-strong galaxies. In Fig. 9 we plot the 'corrected' index given 
by: 

Index = log Mgb + 0.36 log ¥1(3 - 0.07 (8) 

against Vesc and again perform a linear fit to these data. The result- 
ing fit is given by: 

log(Vesc/500kms"^) = 0.16 + 3.571og(Mgb/4A) 

+ 1.291og(H/3/1.6A) (9) 

We find that this fit has a scatter of only a = 0.025 in log Vesc, re- 
ducing the scatter by 22 percent compared to the uncorrected Mgb- 
Vesc relation. This scatter is now consistent with the measurement 
errors. Even with the four galaxies with strongest H/3 removed the 
reduction in scatter is significant, a changes from 0.033 to 0.026, 
a reduction of 20 percent. As a specific example the two galaxies 
at low Vesc lying above the Mgb-Vesc relation fall on the relation 
after this correction is applied. 



4.4 Vesc as an alternative to a and Em 

While a is related to the depth of the potential it is also depen- 
dent on the details of the orbital structure of the galaxy; in galaxies 
with significant rotation or other anisotropy a is a poor tracer of 
whereas the true line-of-sight Vesc is always a reliable measure. 
As an example to support the idea that Vesc is a better predictor of 
local galaxy properties than a we need only look at Emsellem et al. 
(1996). Here the authors study NGC 4594, the Sombrero galaxy, a 
discy edge-on galaxy. In their Fig 21. they show both Mgb vs log a 
and Mgb vs log Vesc- They clearly demonstrate that the local val- 
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Figure 10. Upper figure: Mgb - Vesc profiles for five galaxies from the sam- 
ple. The individual galaxy profiles closely follow the global relation (shown 
as the dotted line). Lower figure: Mgb - a profiles for the same four galax- 
ies. The individual profiles show little resemblance to the global relation 
(again shown as the dotted line), and in some cases exhibit essentially no 
trend with a at all. This illustrates the improvement in using Vesc instead 
of (7 as a tracer of a galaxies dynamical properties. 



ues of Mgb are not significantly correlated with a but are tightly 
correlated with Vesc- 

In our own work we find a similar result. While a is a rea- 
sonable predictor for some galaxies it is generally worse than Vesc, 
and in some cases fails spectacularly to reproduce the line strength 
trends observed with Vesc- This is particularly true of galaxies with 
atypical a maps (central dips in a, counter-rotating discs etc.) In 
Fig. 10 we show the Mgb- Vesc and Mgb-cr relations for a selec- 
tion of galaxies from our sample. As can be clearly seen, while a 
produces a reasonable trend for some galaxies, in others there is no 
observable trend whatsoever yet in these cases the trend with Vesc 
is still quite obvious. We choose to use Vesc because it is a direct 
measure of the potential. 

We also investigated whether Mgb correlates with the local 
surface mass density Em (upper panel. Fig. 11). Hm was calcu- 
lated directly from the MGE models and then scaled using the M/L 
derived from the JAM models. We find that within each galaxy 
there is a tight relationship between Mgb and Em with consistent 
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Figure 11. Upper panel: the Mgb - Em relation for all 48 galaxies. For each 
galaxy the mean Mgb and mean Em of that galaxy has been subtracted 
from each profile. The colour and symbol combination is the same as in 
Fig. 5. Galaxies have the same internal gradients but have very different 
offsets - this relation is local rather than global. The lower panel shows 
the tight relation between Vesc and Em within a galaxy and the lack of a 
connection between different galaxies. 

gradients between galaxies, however, there is no relation between 
the central Em and central Mgb. In this case we find only a lo- 
cal relation; no global relation is apparent. It is interesting to note 
that this implies a constant gradient for Vesc vs Em, but the off- 
set between galaxies shows no such correlation (lower panel. Fig. 
11). This suggests that Em is not physically related to Mgb; the lo- 
cal correlation arising simply because both Em and Mgb decrease 
with galactic radius. Again Vesc appears to be a more significant 
parameter because it shows both a local and global correlation. 

4.5 Influence of a dark matter halo 

As mentioned above we made the assumption that mass follows 
light in order to calculate our Vesc- If the dark matter distribution 
follows that of the stellar mass then this assumption is valid, but 
in general the dark matter profile may be different. The effect of 
a dark halo on the local potential was first explored by Franx & 
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Figure 13. Mgb and Vesc maps for NGC 3377. The solid lines show the ellipses used for extraction of the profiles. Panel a) log Mgb. Panel b) log Vesc derived 
from our best-fitting models. Panel c) log Vesc derived from our thin disc model. As can be seen the isocontours of Vesc produced from the thin disc model 
are significantly flatter than the isophotes or the isocontours of Mgb shown in panel a). 
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Figure 12. Circular velocity (top panel) and escape velocity (bottom panel) 
as a function of radius. The yellow line indicated the total value while the 
red and blue lines represent the contributions due to dark and stellar mass 
respectively. The dashed vertical lines represent the position of 1, 2, 3, 4 
and 5 Re. The change in the gradient of Vesc, between total and stellar 
contribution only, over the region 0-1 Re is 0.07 dex. 



Illingworth (1990). In order to investigate this issue we consider 
the effect of a dark matter halo on a simple galaxy model. We take 
the stellar density to be given by a Hernquist (1990) profile with 
scale radius a = 1 and embed this in a dark halo also represented 
by a Hernquist profile with a = 10: 



p(r) 



M a 



1 



27r r (r + a) 



(10) 



The mass of the dark matter halo was fixed to give a dark mat- 
ter fraction of 50 percent within 5 Re which is consistent with the 
measurements from dynamical studies (see Paper IV, Gerhard et al. 
2001; Thomas et al. 2007) and from lensing (Rusin & Kochanek 
2005; Koopmans et al. 2006). This is similar to the NFW profile 
(Navarro, Frenk & White 1996) in that it has a slope of p ^ 
for r < a, but the Hernquist halo has finite total mass. Recalcu- 
lating our Vesc with this new halo we find that the Vesc gradient 
decreases by 0.07 dex in the interval 0-1 Re (see Fig. 12) between 
the model with a dark halo and that with only a stellar contribu- 
tion. This shows that with these assumptions the dark matter halo 
produces a small but detectable change in Vesc - Over the limited ra- 
dial range covered by our SAURON observations a reasonable halo 
model produces only a modest change of slope. Over a larger ra- 
dial range the dark matter halo can significantly change the slope as 
illustrated for NGC821 and NGC3379 by Weijmans et al. (2009). 



4.6 Mgb and Vesc maps 

While there is a clear correlation between Mgb and Vesc in both a 
global and local sense this is not as obvious when comparing the 
SAURON and Vesc maps. In 18 of the 48 SAURON galaxies the 
isocontours of Mgb are more flattened than those of Vesc (see Paper 
VI. Note, this mostly applies to the Mgb maps. For Fe5015 the 
isocontours are typically rounder than the Mgb contours and so the 
problem is less pronounced if present at all, while for H/3 the maps 
are essentially flat and so not affected by the choice of aperture.) 
Moreover, the Mgb maps show some structure whereas the Vesc 
maps, which trace the potential, are smooth by construction. A lot 
of this difference comes down to rms scatter in the observational 
data which is absent in the modeled Vesc, but the flattening of the 
Mgb isocontours is a significant effect. In particular some galaxies 
(for example NGC3377, see Fig. 13) exhibit a pronounced Mgb 
disc. 
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Figure 14. Edge-on and face-on views of the hyperplane. The directions 
of Vesc and the SSP parameters are shown in the face on view. The red 
symbols on the lower panel show the four galaxies singled-out earlier as 
lying off the Mgb-Vesc relation. The arrows indicate the directions of the 
Vesc and the SSP parameters within the hyperplane. 



Table 3. Principle components analysis 







Age' 
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[a/Fe]' 


Eigen- 
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value 


variance 


PCI 


0.581 


0.286 


-0.214 


-0.731 


2.49 


62 


PC2 


0.531 


-0.388 


-0.606 


0.447 


1.10 


28 


PCS 


0.399 


0.708 


0.275 


0.514 


0.30 


8 


PC4 


0.471 


-0.516 


0.715 


-0.037 


0.10 


2 



Notes: The primed variables are standardised versions of the corresponding 
variables with zero mean and unit variance. The coefficients of the principal 
components are scaled to the variance and sensitive to the range of each 
variable, in the sense that vairables that only vary by a small amount tend 
to have a large coefficient. 



We considered whether the galaxies which exhibit these Mgb 
disc structures would be better fitted by assuming all the Mgb 
comes from a thin disc rather than being uniformly distributed for 
each galaxy. This resulted in only a small change in the determined 
Vesc for these galaxies. The Mgb-Vesc relation derived from the 
disc-based models has a slightly larger scatter than for our best- 
fitting models (a 0.036 for the disc-based values compared to 
cr ^ 0.032 for the best-fitting models) but the relation is still a 
tight one (see Fig. 17). We also considered the effect of extracting 
our Mgb index and Vesc from the maps using a long slit aperture - 
again we recover the tight local and global correlation with similar 
scatter as we found using elliptical apertures. 

While a pure-disc model is clearly unrealistic even this ex- 
treme assumption does not change our main results. We favour a 
scenario of a disc-like structure embedded in a spheroid (see Paper 
VI for further discussion of this idea) to account for the flattened 
Mgb contours and the structure observed in the Mgb maps. Still, 
the key issue here is the link between line strengths and the local 
Vesc, which appears robust against the differing assumptions tested 
above. 

4.7 Vesc and single stellar population (SSP) parameters 

While the line strength- Vesc relations are interesting it is not en- 
tirely clear what they tell us about the formation of local early-type 
galaxies. The measured line strengths are a synthesis of the age, 
metallicity and chemical abundance distributions of the stellar pop- 
ulation. In order to study these more fundamental properties of the 
stellar populations we transform our line strengths into the physi- 
cal parameters, age (t), metallicity ([Z/H]), and alpha enhancement 
([a/Fe]) using the single stellar population models of Schiavon 
(2007). These are not true ages, metallicities and abundances but 
SSP-equivalent parameters assuming each galaxy formed its stars 
in a single burst. While this assumption is clearly unrealistic and we 
should not believe the precise values returned by the model it still 
allows us to make comparisons between the SSP-equivalent values 
for our galaxies. 

The models predict the Lick line strength indices for a wide 
range in age, [Z/H] and [a/Fe] based upon accurate stellar parame- 
ters from library stars and fitting functions describing the response 
of the Lick indices to changes in stellar effective temperature, sur- 
face gravity and iron abundance. The models produce a grid of age, 
[Z/H] and [a/Fe] iso-contours in the Mgb-Fe5015-H/3 space of our 
data. For each data point we find the nearest point on the model 
grids, which gives us the best-fitting age, [Z/H] and [a/Fe] and from 
this we construct age, [Z/H] and [a/Fe] maps. These maps are then 
used to produce the age-, metallicity- and alpha enhancement- Vesc 
profiles using the same method used to construct the line strength 
index- Vesc profiles. This process, along with a more general discus- 
sion of the results is presented in Kuntschner et al. (in preparation). 
Here we confine ourselves to a discussion of the SSP parameters in 
the context of Vesc 

We search for correlations in this four-dimensional Vesc, age, 
[Z/H], [a/Fe] space using principal components analysis (PCA; see 
e.g. Francis & Wills 1999; Faber 1973) the results of which are 
shown in Table 3. As can be seen the first two principal components 
account for 90 per cent of the variance. The properties of local el- 
lipticals are therefore confined to a two-dimensional hyperplane, 
similar to the result found by Trager et al. (2000) but for a instead 
of Vesc Face-on and edge-on views of this hyperplane are shown 
in Fig. 14. There is a lack of points in the bottom right quadrant of 
the face-on view of the plane, due to the upper cut-off in age of 18 
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Figure 15. Upper panels and bottom left: the Vesc vs. [Z/H], Age and [a/Fe]. There is a strong correlation between [Z/H] and Vesc, though not as tight as with 
the individual line strength indices. There is a weaker correlation with age and essentially no correlation with [a/Fe]. None of the figures exhibit the local and 
global correlation observed in the Mgb- and Fe5015-Vesc relations. The sharp cutoff at the top of the upper right and lower left figures is due to the limited 
range of the SSP model. Lower right panel: Edge-on view of the plane connecting [Z/H], age and Vesc, derived from a Hnear fit to the three variables. The 
relationship between the combination of these three variables is much tighter than in the other three panels. Colours and symbols are as described in Fig. 5. 



Gyrs imposed by the SSP model. We checked this result by using 
the models of Thomas, Maraston & Bender (2003) to calculate the 
SSP parameters of our sample and while the precise values in Table 
3 change by 5 percent the conclusion that galaxies are confined 
to a hyperplane is independent of the SSP model used. 

Assuming the hyperplane to be infinitely thin (i.e. the contri- 
butions from PC3 and PC4 are zero) then we can express two of 
our variables in terms of the other two variables. The choice of de- 
pendent and independent variables is entirely arbitrary, but in the 
interests of physical insight we choose [Z/H] and age as our two 
independent variables and seek to express Vesc in terms of them. 
Performing a linear fit to the three-dimensional age, [Z/H], Vesc 
space we find that the variables are related by: 



log 



Ve: 



V SOOkms 



0.85 



+ 0.43 log 



Gyrs 



0.29(11) 



It is important to bear in mind that these are not true ages and metal- 
licities but SSP-equivalent values. This combination of variables is 
shown in the lower panel of Fig. 15. The scatter in this relation is 
greatly reduced from that of any relation between just two of the 



four variables we are considering here as shown in Fig. 15. In Fig. 
16 we show that the local gradients within a galaxy again follow 
the global gradient, though this result is not as tight as the local- 
and-global relation for Mgb- Vesc The global gradient, determined 
from fitting to Re/8 values is 0.79 ± 0.08. The mean of the individ- 
ual gradients is 0.78, with the width of the distribution given by a 
a of 0.40. The typical error on the individual gradients is 0.1 1. The 
global gradient is consistent with the local gradient well within the 
errors. The local gradients again show a broader distribution, im- 
plying an intrinsic scatter of 32 per cent. This is not the case for 
[Z/H] alone; here the local gradients are significantly steeper than 
the global one. The specific combination of age and [Z/H] depends 
on the SSP model used, but the tightness of the plane and the local 
and global connection do not. 
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Figure 16. Local gradients for [Z/H] and 0.85 [Z/H] + 0.43 logt. The histograms show the individual galaxy gradients determined from a linear fit to each 
galaxy's profile. The dotted line shows a distribution with the same mean, a and total area as the individual gradients. The vertical sohd line shows the global 
gradient determined from fitting to Re/8 values only, with the dashed lines indicating the 2 a error. As can be seen the local and global gradients are the same 
for our combination of Z and age but not for Z alone. 



5 DISCUSSION 
5.1 Caveats 

We have already mentioned a few caveats to consider when 
analysing our results. While these points have been discussed more 
fully elsewhere in the text we summarise them here in the interest 
of showing that none of these issues threaten our conclusions. The 
four principal caveats in this work are: triaxiality, inclination, dark 
matter and the shape of the Mgb isophotes. The first three of these 
affect our determination of the Vesc of our galaxies whereas the 
fourth affects the extraction of our Mgb- Vesc profiles. 

i) Bars and triaxiality: Several of our galaxies are triaxial or 
barred systems. While triaxial objects can have very different or- 
bit families to axisymmetric systems the distribution of the matter 
and hence $ will not be significantly different. Furthermore triax- 
ial systems tend to be rounder so the deviations in shape are typi- 
cally small. The M/L is also relatively robust against the assump- 



tion of axisymmetry, which is expected due to the Virial Theorem 
and the tight scaling relations followed by fast and slow rotators, 
e.g. the Fundamental Plane (Djorgovski & Davis 1987). Because 
of this we are able to produce reasonable values for the second mo- 
ments and velocity fields and hence the Vesc of triaxial or barred 
objects under the assumption of axisymmetry. Therefore we do not 
expect that more detailed triaxial modelling (de Lorenzi et al. 2007; 
van den Bosch et al. 2008) will significantly alter our conclusions. 
There is no systematic dependence of the residuals from the Index- 
Vesc relations on bar strength (determined qualitatively by eye for 
our sample) and that our 'clean' axisymmetric sample discussed in 
Section 4 shows no significant improvements in the tightness of the 
relations. 

ii) Dark matter: Our determination of Vesc is based on modelling 
of the photometry and so we are assuming that mass follows light, 
while allowing for a constant dark matter fraction. We note that 
dark matter makes up only a small fraction of the total density in the 
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Figure 17. The Mgb-Vesc relation derived under three different assumptions. Panel a) shows the relation derived assuming all our galaxies are seen with an 
inclination i = 90°. Panel b) shows the results of assuming that all the absorption arises from the equatorial plane of our galaxies, i.e. in a disc. The solid line 
in both panels is taken from Table 2, using the complete sample. As can be seen when comparing this figure with the top panel of Fig. 5 the assumptions made 
in deriving panels a) and b) do not significantly affect our results. 



central regions and hence doesn't significantly affect the potential 
in the region we are studying. There is much observational evidence 
from dynamical studies and from lensing to support this view. In 
Section 4.5 we consider the effect of a dark halo and note that while 
the local gradients in Vesc do change the effect is modest over the 
SAURON field of view. 

iii) Inclination: For many of our galaxies it was possible to esti- 
mate the inclination from methods other than our JAM modelling 
(6 with embedded discs, 16 edge-on objects) and in these cases our 
JAM inclinations fall within the errors on our independent esti- 
mates. For our other galaxies, while we expect that our inclination 
estimates are accurate in most cases (see also Cappellari 2008) we 
also show in Fig. 17 that our Vesc values do not depend strongly on 
inclination. In panel a) of this figure we show the Mgb-Vesc rela- 
tion for our galaxies under the assumption that they are all edge-on. 
As can be seen there is little difference between this panel and the 
top panel of Fig. 5 which shows our Mgb-Vesc relation using our 
best estimates for the inclinations. 

iv) Mgb discs: Finally, as noted in Paper VI, the Mgb isocontours 
do not always follow the isophotes which we are using to extract 
our Mgb-Vesc profiles. We investigated this issue by considering 
the idea that our Mgb absorption comes entirely from a disc and 
re-calculated our Vesc based on this assumption. In panel b) of Fig. 
17 we show the results of this, again, there is little difference be- 
tween panel b) and the upper panel in Fig. 5, though the scatter 
is slightly larger. We also investigated the effect of using differ- 
ent apertures to extract Vesc and Index profiles from the maps by 
varying the ellipticity of our apertures. While the scatter increased 
slightly when circular apertures were used the relations did not sig- 
nificantly change. 

5.2 The line strength-Vesc relations 

It is clear from the tightness of the correlations shown in Fig. 5 that 
the stellar populations of early-type galaxies are closely linked with 
the depth of the local potential they reside in (characterised in this 
study by Vesc). That this should be the case is by no means obvious. 
While we might expect that the formation of a star is influenced by 
the potential that it forms in it is perfectly possible for that potential 



to have changed significantly between the star's formation and the 
present day. We expect that the availability of gas and it's ability 
to cool will also play a role. The tight correlation observed is more 
easily accommodated in a monolithic collapse scenario for star for- 
mation in which the gravitational potential $ is largely unchanged, 
but it is clear that galaxies do not form in this way - we need to 
address how mergers fit into this picture. 

While in the monolithic collapse scenario the potential does 
not change the same is not true of mergers; in this case the potential 
can be significantly altered as more mass is added to the galaxy and 
the distribution of that mass can also be changed. In this case it is 
clear that the potential a star forms in and the potential we observe 
it in several billions years later are different. It has been known for 
some time (White 1978; Barnes 1988) that during a merger the stars 
are preserved in their rank-order of binding energy, in the sense that 
the most deeply bound stars before the merger are also most deeply 
bound after the merger. This suggests a possible link between the 
potential a star formed in and the potential it finds itself in after 
the merger. More recent work by Hopkins et al. (2008) has argued 
that in both wet and dry mergers the radial gradients of metallicity 
are preserved, again suggesting that the tightness of the observed 
line strength-Vesc relationships can be compatible with hierarchical 
merging. However it is the detail that the local and global relations 
are the same that is our key result and more detailed modelling is 
required before we can properly compare model predictions to our 
result. 

Four galaxies in our sample show signs of recent star forma- 
tion, and these galaxies are significant outliers in the Mgb- and H/3- 
Vesc relations. But, when we convert our line strength indices into 
SSP-equivalent parameters we find that these four galaxies are con- 
fined to the same region in the four-dimensional parameter space 
of Vesc, t, [Z/H] and [a/Fe] as those galaxies which lie on the 
Index- Vesc relations. As these objects all have unusually high H/3 
they are all likely to be relatively young objects. This suggests that 
even recently disrupted objects have some regular properties con- 
nected with the potential $ that survive whatever process moved 
the galaxy off the Index- Vesc relation. It seems likely that as these 
objects age they will return to the Index- Vesc relations. 

It is interesting to ask why we observe a tight local and global 
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Figure 18. The line strength - Vesc relations predicted by re-inverting the 
SSP model of Schiavon (2007) and applying the constraint of Equation 1 1 
at constant [a/Fe/] = 0.33 and minimum age, t > 2 Gyrs. The grid of points 
show the predictions of the SSP model. The red and blue lines show lines of 
constant age and [Z/H] respectively. The green line is the best fit to the ob- 
servations, given in Table 2. In the lower panel the grey shaded area shows 
the region occupied by the observations. With these constraints the SSP 
model is able to reproduce the observed Mgb and H^^-Vesc relations but the 
condition of fixed [a/Fe] is required to reproduce the Fe5015-Vesc relation. 



relation in the Mgb-Vesc relation but not in the case of the other two 
indices. In order to investigate this issue we begin with our SSP hy- 
perplane equation, Equation 11. We re-invert the SSP model grid of 
Schiavon (2007) and consider the constraints on the line strength- 
Vesc relations implied by Equation 1 1 , using a constant [a/Fe] of 
0.33, the mean value for galaxies in our sample. We also impose 
a minimum age of t > 2 Gyrs. This is shown in Fig. 18. By de- 
sign the Mgb-Vesc relation produced in this way tightly follows the 
observations. In the case of H/3 the region allowed by the model 
(the grid of points) is well matched to the region occupied by the 
observations (shaded area). The observed local gradients broadly 
follow the lines of constant age in the modeled region, which may 
suggest why the local and global connection is not observed be- 
tween H/3 and Vesc In the case of Fe5015 the model predicts a 
tight correlation with Vesc but one that is somewhat steeper than 
the observed relation. Allowing other values of [a/Fe] is required 
to match the observations; this is consistent with the SSP hyper- 
plane having a small thickness. We again note that the steeper local 
gradients in Fe5015 appear to follow lines of constant age. These 
results suggest two possible conclusions; either the local and global 
correlation in Mgb is a conspiracy of the interaction between age 
and metallicity in producing stellar absorption indices or that there 
still remain weaknesses in the SSP models used to derive Equation 
11. 

We find essentially no difference between the relations for 
fast- and slow-rotators, which are thought to have significantly dif- 
ferent assembly histories. There are clear differences between the 
two classes of galaxies in many of their properties (see Emsellem et 
al. 2007, for a discussion of these differences) but not in the Index- 
Vesc relation. What does this tell us about the assembly of these 
galaxies? If fast rotators are the progenitors of slow rotators and 
we assume that fast rotators lie naturally upon the Index- Vesc rela- 
tions we have observed whatever process leads to these differences 
between fast- and slow-rotators must preserve the links between 
stellar population properties and the gravitational potential. If, as 
Di Matteo et al. (2009) suggest, equal mass dry mergers between 
giant elliptical galaxies can significantly alter the metallicity gradi- 
ents of the remnant, the lack of a difference between the relations 
for fast- and slow-rotators may provide a significant constraint on 
the modelling of the formation of these galaxies. 



6 CONCLUSIONS 

In this work we have examined the link between the local escape 
velocity, Vesc (determined from photometric observations and dy- 
namical modelling) and the local line strength indices. We discuss 
the impact on our results of non-axisymmetry, dark matter, inclina- 
tion and substructure within the line strength maps. Single stellar 
population models were used to convert our line strength measure- 
ments into representative values for the age t, metallicity [Z/H], 
and alpha enhancement [a/Fe]. We then used some simple models 
to explore the impact of the observed correlations on the formation 
history of early-type galaxies. The main findings of this work are 
as follows: 

i) The line strength indices Mgb and Fe5015 are correlated with 
Vesc (both with rms scatters of 0.033) while H^^ is anti-correlated 
with Vesc (with an rms scatter of 0.049). Using the models of Schi- 
avon (2007) the scatter in the Mgb relation corresponds to a spread 
A[Z/H] 0.2 at a fixed age of 9 Gyrs. The scatter in the H/3 re- 
lation corresponds to a spread At ^ 4.7 Gyrs with [Z/H] fixed at 
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solar metallicity. The tightness of these relations provide an impor- 
tant check for simulations of early-type galaxy formation. (In com- 
parison the index - cTe relations have rms scatters of 0.028, 0.030 
and 0.046 for Mgb, Fe5015 and H/3). 

ii) For Mgb the correlation within a galaxy (the local relation) is 
the same as that between the central values of different galaxies 
(the global relation). This is the key difference when considering 
Vesc compared to using a. 

iii) For outhers characterised by high H/3 the residuals in the Mgb- 
Vesc relation correlate with H/3. We use this correlation to mod- 
ify our Index- Vesc relation such that log(Vesc/500kms~^) = 
0.16 + 3.571og(Mgb/4A) + 1.29 log (H/3/ 1.6 A). The scatter of 
this corrected relation is consistent with the measurement errors. 

iv) We divided our sample into several sub-populations: SOs and 
ellipticals, field and group/cluster objects and fast- and slow- 
rotators. The Index- Vesc relations for each of these sub-populations 
are consistent with the relations for the entire sample. We find no 
dependence on these simple divisions into morphological type and 
environment. We also find no significant difference between fast- 
and slow-rotators. 

v) When converting these line strength measurements to SSP pa- 
rameters we find that all the galaxies are confined to a two- 
dimensional plane within the four-dimensional space of Vesc, 
age [Z/H] and [a/Fe]. This plane is described by the equation: 
log(V;sc/500kms-^) = 0.85[Z/H] + 0.43 log(t/Gyrs) - 0.29. 
Those galaxies that were outliers in the Index-Vesc relations do not 
stand out in this SSP-hyperplane. 

vi) We find that in the Z-Vesc diagram the local gradients are sig- 
nificantly steeper than the global relation. When we consider the 
above combination of Z and age we recover the local and global 
relation, in that the local gradients are the same as the global one. 
This tight relation does not depend on the SSP model used. 

How the connection between stellar populations and the gravi- 
tational potential, both locally and globally, is preserved as galaxies 
assemble hierarchically presents a major challenge to models. 
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and red symbols indicate slow rotators. The solid line is a fit to the data, 
excluding the galaxy with the lowest M/L, NGC3489. 
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APPENDIX A: MASS -TO-LIGHT RATIO 
CORRELATIONS 

In Paper IV we examined the correlation between the dynamical 
mass-to-light ratio (M/L)jeans, and several other properties of lo- 
cal early-type galaxies, including erg, H/3 and the virial mass-to- 
light ratio, (M/L) Virial- Here we revisit those correlations with the 
additional 24 galaxies contributed by this work. The analysis was 
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Figure A2. (M/L) jeans vs (M/L)viriai • The solid line is a fit to the data 
while the dotted line shows a one-to-one correlation. Symbols as in Fig. 
Al. 
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Figure A3. Dynamical M/L versus the observed line strength index, H/3. 
Symbols as in Fig. Al. The shaded region indicates the predictions of the 
SSP models of Vazdekis et al. (1996) and Vazdekis (1999) using a Saltpeter 
IMF 



carried out precisely as described in Section 4.2 of Paper IV, us- 
ing mass-to-light ratios derived from JAM modelling rather than 
Schwarzschild modelling. Where available we used surface bright- 
ness fluctuation distances, otherwise distances are redshift only and 
have a significantly larger error (Tonry et al. 2001; Paturel et al. 
2003; Mei et al. 2007). All mass-to-light ratios were converted to 
/-band using galaxy colours obtained from Tonry et al. (2001) and 
Prugniel & Heraudeau (1998). As in Paper IV we adopted an error 
in cTe of 5 percent and a 6 percent modelling error in (M/L)jeans 
to which we quadratically co-added the distance errors. The fit 
was carried out by quadratically adding an intrinsic error to make 
X^/z^ — 1, where ly is the degrees of freedom. To preempt our con- 
clusions we find no significant change from the results presented in 
Paper IV. 

We find a tight correlation with (Je, with an observed rms scat- 



ter of 30 percent, shown in Fig. Al. This implies an intrinsic scatter 
of 1 1 percent. The best fitting relation has the form: 

/ ^ \0.89±0.09 

(M/L)._ = (3.77 ± 0.14) (^^^^^) (Al) 

We also compared our dynamical (M/L) jeans to the observed virial 
(M/L)viriai = f3Recre/GL (whcrc we adopt /3 = 5 as found in 
Paper IV), shown in Fig. A2. With our adopted modelling error 
of 6 percent in the (M/L) jeans the scatter in (M/L)viriai required 
to make x^/^ — 1 is 26 percent. The best fitting relation has the 
form: 

(M/L)jeans « (M/L)^fif (A2) 

We also present the relation with H^^ (see Fig. A). This is closely 
related to the population mass-to-light ratio (M/L)pop, which is 
largely driven by variations in H/3. We again reproduce the trend 
found in Paper IV. We also indicate the predictions of the SSP mod- 
els of Vazdekis et al. (1996) and Vazdekis (1999) as shown in fig. 
16 of Paper IV 



APPENDIX B: PSF PARAMETERS USED IN THE MODEL 
MGES 

Table Bl contains the details of the point spread functions used 
when evaluating our model MGEs. These PSFs were calculated by 
fitting circular Gaussians to PSFs produced using the Tiny Tim soft- 
ware. The PSFs have the form: 

PSF = E^^iGfcexp[-i?V2a^]/2^a^ 



APPENDIX C: FURTHER EXAMPLES OF VELOCITY 
FIELDS CALCULATED FROM JAM MODELLING 

In Fig. C we show the results of JAM modelling for eight galaxies 
from our sample. Further examples are shown in Cappellari (2008). 
Details of the JAM modelling can be found in Section 3.2. We show 
the observed and modeled results for both the first and second mo- 
ments of the velocity for each galaxy. As can be seen in all cases we 
achieve good agreement between the observations and our models. 



APPENDIX D: MULTI-GAUSSIAN EXPANSION 
PARAMETERS FOR 24 GALAXIES IN OUR SAMPLE 

In Appendix B of Paper IV we presented the (distant-independent) 
MGE parameters for 23 galaxies from the SAURON sample^. In 
Table Dl and D2 we present the MGE parameters for the remaining 
23 galaxies. The constant-PA models were obtained by fitting the 
HST photometry (where available) at small radii and the ground- 
based MDM photometry at large radii. They provide an accurate 
description of the surface brightness of the galaxies from ^ 0.01 
arcsec to about twice the maximum aj used in each galaxy (usu- 
ally corresponding to 5 — 10 Re). Dust and bright foreground stars 
were excluded from the fit. The matching of the different images is 
described in Section 3.1. The deconvolved surface brightness E is 
defined as follows: 

^ The MGE parameters for NGC 2974 were given in Krajnovic et al. 
(2005). Those for NGC 3032 were given in Young, Bureau & Cappellari 
(2008). 
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Figure CI. Examples of the first and second moments of the velocity fields, v (/ii) and a/u^^F^ (/i2) derived from our JAM models alongside the bi- 
symmetrised measured /ii and 112 from the SAURON kinematics. From top to bottom for each galaxy the figures show: observed /i2, model 112, observed /ii, 
model /ii . We show here a selection of galaxies from our sample to illustrate the range in the quality of fits we achieved. 
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Table Bl. Parameters of the model MGE PSFs 



Instrument/Filter 
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Notes: The model PSFs are formed from the sum of circular 
Gaussians fitted to TinyTim PSFs and have the form: PSF = 
Tl^^^Gk exp[— i^^/2(7^]/27r(7^. The numerical weights are normalised 
such that S^^-L 1 

where the model is composed of N Gaussian components of dis- 
persion (jj, axial ratio qj and peak intensity Ij. The coordinates 
(x^ y^) are measured on the sky plane, with the x-SLxis correspond- 
ing to the galaxy major axis. The total luminosity of each Gaussian 
component is given by Lj = 27vljajqj. See Cappellari (2002) for 
details. The obscured areas in each figure were masked in the fitting 
process. 
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Table Dl. MGE parameters for the deconvolved /-band surface brightness 
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-1.422 


0.910 


5.049 


-1.589 


0.692 


4.148 


-0.728 


0.770 
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10 
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0.636 
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0.700 


1.818 


2.154 


0.780 



Note: * indicates where a negative Gaussian was used to achieve a satisfactory fit. This only occurs in the most discy objects. 



ne SAURON Project - XIV 25 



Table D2. MGE parameters for the deconvolved /-band surface brightness 
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Note: * indicates where a negative Gaussian was used to achieve a satisfactory fit. This only occurs in the most discy objects. 



